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< ^BSTRACT: We study the injection of air in a granular suspension. We use a linear Hele-Saw cell filled with 
a suspension which is displaced by air, leading to a Saffman-Taylor (fingering) instability. For the suspension, 

^ve use an iso-dense mixture where the fluid and the particles have the same density. The volume fraction of 
particles can thus be adjusted over a wide range. We discuss the question of an effective rheology inside the cell 
well as the pattern formation as a function of the granular compacity. We finally report results on the finger 

~gvidth for stable fingers and the thresholds for their destabilization. 
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Q] INTRODUCTION 

^Viscous fingering has received much attention as an 
f^rchetype of pattern-forming systems and was stud- 
Qed intensively for Newtonian fluids (Saffman & Tay- 
1958, Bensimon et al. 1986, Homsy 1 987) and 
non-Newtonian fluids (L indner et al. 2002b . It is an 
^mportant model system suited for a better under- 
standing of fluid injection in porous media with im- 
CXportant practical applications such as crude oil recov- 

q ■ Here, we consider a case where the fluid to be dis- 
iTtolaced is an isodense non-Brownian suspension. The 
^volume fraction of the grains can be easily controlled 
^nd thus, the rheology of the suspension can be var- 
Cjed over a wide range: going from an effective vis- 
icous fluid up to a jammed dense paste. This system 
r W;an be seen as a fundamental tool to understand many 
g;omplex injection problems such as those existing in 
weakly consolidated porous media. Its potential inter- 
^sest would lie in a better understanding of many prac- 
tical processes such as the stabilization of saturated 
^oils by air injection or the consolidation of oil ex- 
traction wells. 

A difficult question is still to understand the interfa- 
cial dynamics between a suspension and another fluid. 
Viscous fingering is mainly governed by the local vis- 
cosity at the finger tip and is also very sensitive to the 
presence of particles. We thus expect that studying the 
pattern formation will lead to a better understanding 
of particles dynamic at the interface between suspen- 
sion and air (or another injection fluid). 

The paper is organized as follows. In Sec. 2 we will 
recall the basic equations for the Saffman-Taylor in- 



stability. Sec. 3 describes the set-up and experimental 
methods. In Sec. 4, the experimental results concern- 
ing the rheology of suspensions and the Darcy's law 
are presented and discussed. In Sec. 5, we interest our- 
selves at the finger shape as well as the finger width. 
Sec. 6 gives a summary of the obtained results 

2 THE SAFFMAN-TAYLOR INSTABILITY 
2.1 Presentation of the instability 
The Saffman-Taylor instability is typically studied in 
a thin linear channel or Hele-Shaw cell (Fig. [II). The 
width of the cell W is chosen to be large compared to 
the channel thickness b. The cell is filled with a vis- 
cous fluid which is then pushed by air. The properties 
of the viscous fluid are its viscosity /i, its surface ten- 
sion 7 and its density p. The viscosity and the density 
of air are neglected. 

When air pushes the viscous fluid (for example due 
to an imposed pressure gradient VP), the interface 
between the two fluids destabilizes. This destabiliza- 
tion leads to the formation of a viscous finger of width 
w propagating at velocity U. In general, one considers 
the relative finger width: A = w/W. 
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Figure 1 : Schematic drawing of the experimental set-up. 
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2 . 2 Averaged Darcy 's law 

For Newtonian fluids, the motion in the Hele-Shaw 
cell is described by the two-dimensional gap averaged 
velocity field u. It is given by the local Darcy's law, 
which relates the local pressure gradient to the fluid 
velocity. 

Away from the interface the flow can be considered 
as uniform and one can relate the mean flow velocity 
V to the pressure gradient VP: 



V 



12r7 



-VP. 
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2.3 Finger selection 

The relative finger width A is determined by the cap- 
illary number Ca = rjU/j which represents the ra- 
tio of viscous forces and capillary forces. The viscous 
forces tend to narrow the finger, whereas the capillary 
forces tend to widen it. One thus anticipates that the 
relative width of the viscous fingers decreases with 
increasing finger velocity. For large values of the cap- 
illary parameter, A reaches a limiting value of about 
half the channel width. 

The control parameter of the instability is 1/5 = 
\2{W/b) 2 Ca with W/b the aspect ratio of the Hele- 
Shaw cell. When scaled on 1/B, measurements of A 
for different systems all fall on the same universal 
curve (Saffman & Taylor 1958, McLean & Saffman 
1981). 

3 EXPERIMENTAL SET-UP 
We use polystyrene spheres of 40 /im in diame- 
ter in a modified silicone oil DC704 of viscosity 
r]o = 59.5 mPa.s and surface tension 7 = 31.0 mN/m. 
The density p of the fluid matches the one of the 
grains (1.05 g/cm 3 ) closely leading to what we call 
an iso-dense suspension. This allows us to control the 
fraction of grains over a wide range and no sedimenta- 
tion occurs on the timescale of our experiments. Here 
we work with grain fractions (i.e. the amount of grains 
on the total amount of suspension) $ up to 30 %. 

We worked in a linear Hele-Shaw cell consisting 
of two glass plates separated by a thin Mylar spacer 
(Fig. [l]). These plates were horizontal and clamped 
together in order to obtain a regular thickness of the 
channel b of 0.75 mm for a width W of 4 cm. The cell 
was filled with the suspension and compressed air was 
used as the less viscous driving fluid. 

Fingers or more complex patterns were captured 
by a CCD camera, coupled to a data acquisition card 
(National Instruments) and a computer. In particular, 
this allowed for measurements of the relative width 
A = w/W as a function of the velocity U for stable 
fingers. For each grain fraction several experimental 
runs (between 20 and 30) were performed using dif- 
ferent applied pressure drops leading to many differ- 



ent finger velocities. 

The rheology of our suspensions was studied in a 
double Couette geometry (2x 0.5 mm) on a Thermo- 
Haake rheometer. 

4 RESULTS : RHEOLOGY AND DARCY' S LAW 
4.1 Rheology 

Figure |2] shows the viscosity of the suspensions as a 
function of the shear rate 7 for different grain frac- 
tions $. 

For all suspensions one observes Newtonian behav- 
ior over a wide range of shear rates (ls _1 <7<100 
s _1 ). Our fingering experiments are indeed performed 
in this range of shear rate and we can thus define a 
constant viscosity rjBh from the rheological measure- 
ments for each suspension. 

For large shear rates (7 > 100 s _1 ) and high grain 
fraction ($ = 30 %) weak shear thinning is observed. 

For low shear rates (7 < 1 s _1 ) and high grain frac- 
tions a decrease of the viscosity is observed which 
could be due to the existence of a yield stress, whereas 
for low grain fractions an increase of viscosity is ob- 
served. 

Note that these results are in qualitative 
and quantitative agreement with previous work 



(Zarrag a et al. 20 00). 
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Figure 2: Rheological measurements: viscosity as a func- 
tion of the shear rate. 



4.2 Darcy's law 

Assuming the flow far away from the finger to be uni- 
form, we expect that the classical Darcy's law (Equa- 
tion 1) linking the gap averaged fluid velocity V to the 
imposed pressure gradient VP in our Hele-Shaw cell 
remains valid. 

The imposed pressure gradient is calculated by: 
VP = Ap/L where Ap is the applied pressure drop 
and L the distance between the finger tip and the exit 
of the cell. Mass conservation allows to obtain the 
velocity V of the fluid far away from the interface 
from the finger velocity U simply by using V = XU 
(Lindner et al. 20001) . 
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Following equation 1, results are described on Fig- 
ure ^representing V.(l2r) Rh /b 2 ) as a function of VP. 
If one sompares with what is obtained from Darcy's 
law (Fig. |3 one observes systematic deviations to- 
wards higher slope values, which increase with in- 
creasing volume fraction. This might indicate that 
there exists an effective viscosity rjc in the cell which 
is lower than rj Rh obtained with the rheometer. 
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Figure 3: Darcy's law for the suspensions. The straight 
line represents Darcy's law expected in a Hele-Shaw cell. 



4.3 Flow in the cell 

Figure |4]reports the relative difference of the viscosi 
ties: 



-($) 



(2) 



VRh ' VRh($) 

as a function of $. The cell viscosity r] C is always 
lower than r] Rh and the relative difference increases 
with the grain fraction. 

These results could be qualitatively explained by 
a granular structuring of the flow in the thickness of 
the cell fLyon & Leal 19 98): on one side, the flow is 
no longer parabolic but evolves towards a plug flow 
with increasing grain fraction and, on the other side, 
the grain concentration is no longer homogeneous and 
becomes higher in the central region at low shear. It is 
clear that such a restructuring of the flow could lead to 
a lower average viscosity rjc- Unfortunately, no direct 
link with an effective cell viscosity has been estab- 
lished so far. 
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Figure 4: Relative viscosity difference as a function of 
grain concentration <3?. 



5 RESULTS: FINGERS AND PATTERNS 

5 . 1 Experimental observations 

Figures |5k,b,c,d,e,f present the typical evolution of 
patterns observed for increasing velocity, valid for 
most grain fractions ($ up to 40 %, not presented 
here). 
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Figure 5: Typical evolution of the finger shape with in- 
creasing velocities: 0.3 (a), 2 (b), 5 (c), 10 (d), 12 (e) and 
22 cm/s (f). 



For low velocities, the finger is stable and symmet- 
rical with respect to the cell central axis (Fig.|5k). For 
higher velocity, it become asymmetric (Fig. |5b), al- 
though it remains stable (i.e. we could measure its 
width). 

The first destabilization occurs on one side of the 
finger like a sinusoidal wave (Fig.|5t). Thereafter, for 
increasing velocity, this type of destabilization occurs 
on both sides (Fig. EJl). 

At ever larger velocity, we observe classical 
destabilization of Saffman-Taylor fingers like side- 
branching (Fig.|5k) or tip-splitting (Fig. Or). 

5.2 Fingerwidth 

For stable fingers (Figs [5^, b), we measured the rela- 
tive finger width A. When representing A as a function 
of the control parameter l/B of the instability one has 
to decide which viscosity to use. 

Figures |6] and |3 represent A as a function of l/B 
calculated with r] Rh (Fig . (6J) and with i] C (Fig. 13). The 
line represents a fit to the data for the pure fluid and 
gives thus the master curve of the instability. 

We observe an increasing dispersion of the results 
with grain fraction. However, on Figure the data is 
scattered around the master curve whereas on Figure 
13 all data seem to be below this curve. 

As the finger width depends on the viscosity at the 
finger tip, one expects that a non-uniform flow and 
grain distribution will affect the finger width. Thus, 
one could interpret these results in terms of an ef- 
fective viscosity at the finger tip which would vary 
around r] Rh with rjc being the lower limit. 
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Figure 6: Finger width A as a function of the control pa- 
rameter 1/B using Tjjih. The straight line represents a fit to 
the data for the pure fluid. 
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Figure 7: Finger width A as a function of the control pa- 
rameter 1/B using tjq. The straight line represents a fit to 
the data for the pure fluid. 



5.3 Instability thresholds 

Following the evolution of the finger shapes (Figs 13), 
one observes several instability thresholds. On Fig- 
ure HI the first threshold (apparition of oscillations on 
one side Fig. |3t) and the threshold for side branch- 
ing (Fig.|5k) are plotted. Here r]c is used to obtain the 
value of 1/B C . 

We observe little variations of the first threshold 
with $ which is found to be 1/B C « 2500 whereas 
in the case of pure fluid (i.e. $ = 0), we observe stable 
fingers up to 1/B « 10,000. 

Following the work of Bensimon et al. (1986), 
one can obtain a critical amplitude of the perturba- 
tions leading to a finger destabilization at a value 
1/B C = 2500 for a finite amplitude A c » 5~20 fim. 
Furthermore, this work predicts the destabilization to 
be independent on the wavelength of the perturbation. 
These predictions are in good agreement with our ex- 
perimental data if we consider the size of the grains 
(40 /mi) to be of the order of A c and the threshold to 
be weakly dependent of fraction of the grains (i.e. the 
concentration of grains at the finger tip) which should 
be proportional to the wavelength of the perturbation. 
To test this further experiments are planned changing 
the particle diameter. 
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Figure 8: Instability thresholds: □, apparition of oscilla- 
tions on one side of the finger; •, apparition of side branch- 
ing. 



6 CONCLUSION AND PERSPECTIVES 
In summary we have reported results on the injec- 
tion of air in granular suspensions. The dynamics of 
the air/suspension interface can be described by a 
Saff man-Taylor like instability. Rheological measure- 
ments and the use of Darcy's law for flow in the Hele- 
Shaw cell lead to different values for the viscosities 
of the suspensions. This could be explained by a non- 
homogeneous flow in the cell. 

Furthermore we report the evolution of the ob- 
served patterns with finger velocity. For stable fingers, 
an increasing scatter of the data of the relative fin- 
ger width with the grain fraction was observed, which 
might be linked to the local viscosity at the finger tip. 

Finally, an analysis of the instability thresholds 
leads us to underline the importance of the grain size 
for the destabilization of the fingers. 

In further studies of the dynamics of the interface 
between a pure fluid and a suspension, two cases 
seem to be of interest: the case of dense suspensions 
(> 50 %) and the case of miscible flow where the 
pushing pure fluid and the suspending one are iden- 
tical. 

We thank M. Cloitre (ESPCI, Paris) for helping us 
with the rheological measurements and Daniel Bonn 
(ENS, Paris) for helping us with the experimental set- 
up. 
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